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We present thermodynamic measurements of various phys- 
ical observables of the two dimensional S=l/2 isotropic quan- 
tum Heisenberg antiferromagnet on a square lattice, obtained 
by quantum Monte Carlo. In particular we have been able 
to measure the infinite volume limit of the uniform suscep- 
tibility up to the inverse temperature f3 = 40, the analysis 
of which reveals the correct asymptotic behavior in excellent 
agreement with the prediction of chiral perturbation theory. 
The issue of the existence of a crossover from quantum critical 
to renormalized classical regime is clarified. 
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The square lattice quantum Heisenberg antiferromag- 
net (QHA) is important in condensed matter physics 
since it describes the critical behavior of undoped insu- 
lating parent compounds of the high T c superconductors. 
Some materials believed to be well described by these 
models are La 2 Cu0 4 0, Sr 2 Cu0 2 Cl @, La 2 Ni0 4 f§ and 
K 2 NiF4 ||. Most existing theoretical treatments of the 
long wavelength, low energy behavior of the QHA are ac- 
tually based on an effective field theory, the (2+l)D 0(3) 
non- linear a model (NLctM) §]-§. 

Although the mapping from the QHA to the NLoM 
may be justified on the general ground of universality, it 
is rigorous only for sufficiently large values of the mag- 
nitude of the quantum spin (S). There is also a subtle 
problem due to the existence of Berry phase terms, which 
are present in the effective field theory of the QHA 0, 
but not in the NLcrM. This term is believed not to be 
relevant |]|| , but this is still a matter of debate |]] . 

The results of previous experimental and numerical 
studies on the S=l/2 QHA show a leading exponen- 
tial temperature dependence of the correlation length, 
in agreement with the prediction of the theories [^^). 
There nevertheless exist some systematic discrepancies 
between the observed values and theoretical predictions: 
The two-loop order formula of the correlation length de- 
viates by more than 15% from both measurements [^|J^] 
and numerical simulations |[(J[ll). The deviation be- 
comes worse as the value of S increases jllj, in contrast 
to naive expectations. Recent neutron scattering exper- 
iments on SriCuOiCli M, which is regarded as a more 



from the theoretical prediction for the peak value of the 
static structure factor. The presence of a crossover from 
quantum critical to renormalized classical behavior pre- 
dicted by the theories remains controversial with 
differing claims reported in the literature ]lO| , p^ ,p[ . 

In this Letter we present results of a very large scale 
quantum Monte Carlo study on the S=l/2 square lattice 
QHA up to the linear size of the lattice L=1000. Our 
study resolves the above discrepancies, and clearly con- 
firms the validity of the theory . Using the continu- 
ous time version |l^] of the loop algorithm jl4j that elim- 
inates the systematic error due to finite Suzuki- Trotter 
number, we measure the thermodynamic values (infinite 
volume limit values) of various physical observables such 
as the uniform susceptibility (x u ), staggered suscepti- 
bility (Xs)i second moment correlation length (£), peak 
value of the staggered structure factor Sq at Q = (it, n) 
and internal energy (£). For each T and L we performed 
10 6 ~ 10 7 single loop updates after thermalization. 

In order to monitor finite size effects in our measure- 
ments, we repeated the measurements for each temper- 
ature with varying lattice size and found that the mea- 
sured values of £, \ s and Sq become size independent 
under the condition L/£ > 7, within the typical relative 
statistical errors of 0.3 percent or better. This condition 
is very restrictive for the measurements of the infinite 
volume value of those quantities. It turns out however 
that Xu and £ can be measured reliably on much smaller 
lattices. The uniform susceptibility could be measured 
up to /3 ee J/T = 40, where \ is of order O(10 19 - 10 20 ). 
Varying L from 20 to 120 we found that the data are 
already size independent for L > 80 within the typical 
statistical error of 0.3 percent or better. 

A selection of our data is shown in Table 1. The com- 
plete data will be presented in a forthcoming publication 
p"5[ . In the following we discuss our measurements and 
fit them to theoretical predictions. 

Uniform susceptibility: Chiral perturbation theory ||, 
for Xu as T — > 0, predicts (with convention J — Ti — 
cj[iB = = 1 hereafter) 
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ideal realization of the S=l/2 QHA than the previously 
used compound La^CuOi 111], show systematic deviations 
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where x± = Ps/c 2 - It turns out that the fit of our data to 
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FIG. 1. The uniform susceptibility \u versus T. The inset 
shows the low-T region enlarged. The fit to the low-T predic- 
tion Eq. ([!]) is drawn as a solid line, the fit to the quantum 
critical linear behavior with is shown as a dashed line. In the 
latter fit the slope is completely determined by the value of 
the Xu(0)i extracted from the first fit. 

Eq. (|l|) becomes stable only for data with f3 > 4.5, with 
X 2 /Ndf — 0.6. The estimated values of the parameters 
from the fit arc 



p s = 0.179(2), c= 1.652(1). 



(2) 



These values are in remarkably good agreement with 
those obtained from spin wave theory of the QHA []l6| . 
Our estimate of \± is 



Xx = 0.0656(1), 



(3) 



which is also consistent with the result of the spin wave 
theory @. 

Also predicted theoretically || is a crossover from 
quantum critical behavior to renormalized classical one, 
that is, from [Bl 



Xu = ^{A U T + B u (p s )\ 



(4) 



to Eq.(|l]). In Eq. (|J) the constant A u is universal. The 
best estimate from QMC is A u = 0.26 ± 0.01 g. B u is 
only known to leading order in a 1/iV expansion Q as 
B u k. 0.57p s . We observe (Fig. 1) that in a reasonably 
broad range 0.3 < T < 0.5 Xu is linear in T with the 
expected slope A u and an offset B u « 0A7p s , reasonably 
close to the analytical estimate. 

Internal Energy: We find that the energy £ becomes 
size independent under the condition L/£ > 3. Our mea- 
surements are over 0.25 < (3 < 5.5, usually by varying f3 



FIG. 2. Ratio of the measured and predicted correlation 
length, i.e., ^/^hn(Ps, c) with the values of p s and c given in 
Eq. (2), as a function of T. 



by 0.25, on lattices of size up to L=1000. The theoretical 
prediction given in Ref. || is 



£(T) = E a + E 3 T 3 



E,T\ 
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Due to considerable uncertainties in three parameter fits, 
we here fix the value of E% with the value of c given from 
the fit of Xu, i-e, c = 1.652 rather than treating E% as a fit- 
ting parameter. It turns out that only data for (3 > 4.25 
fit reasonably well to Eq. (||) with x 2 /%f — 1-0. We 
observe however that the fit is still slightly unstable in 
the sense that the values of the fitting parameters change 
mildly with the range of T selected for the fit. The 
ground state energy we have extracted is 



E = -0.66953(4), 



(0) 



which may be compared with Eq = —0.6693 ~ —0.6694 
obtained from the ground state properties of QHA 
[|Tt| , P~8|| . Due to the mild increasing tendency of E with 
restricting the fit to lower temperatures, our estimate 
(Eq. (||)) may be regarded as a weak lower bound of the 
correct value. 

Correlation length: The prediction of £ from the two 
loop order of the chiral perturbation theory up to the 
first order of T reads lal 
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We measured the infinite volume limit £ over the in- 
verse temperature range 0.25 < (3 < 4.95, corresponding 
to 0.289(2) < £ < 120.5(4), and fitted our data to Eq. 
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(0) in various temperature ranges . We find that the fits 
are unstable, i.e., the value of p s (c) decreases (increases) 
systematically as the data at higher T are removed in the 
fit. Considering data with £, m i n — 39.2(1) in the fit we 
find 



0.185(1), c= 1.442(3). 



(8) 



Due to the systematic tendency this value of p s (c) should 
be regarded as the upper (lower) bound of the correct 
asymptotic value. Although those extracted from our \u 
data , Eq. (§), are indeed consistent with the the bounds, 
our £ over 0.25 < (5 < 4.95 strongly deviate from the 
asymptotic expression ( FigJ^). However we wish to note 
that both our data and the experimental measurements 
(H can be fitted quite well by Eq. (0) with the correct 
p s if one leaves the prcfactor a free fitting parameter, as 
was done in Ref M . 

The theory ^,^also predicts a crossover from the be- 
havior in the quantum critical regime, given by c/£(T) = 
AqcPs+BqcT, to the asymptotic behavior Eq. (0). The 
expressions for Aqc and Bqc were obtained to one-loop 
order of the renormalization group calculation A 
previous Monte Carlo study |l(J claimed that all £ fit a 
simple exponential form even for (3 as low as 0.25, con- 
trary to our data. They thus concluded that there is no 
crossover. On the other hand a series expansion study 
for £ < 10 claimed to see a crossover We indeed ob- 
serve that l/£ is linear in T over 1.273(6) < £ < 3.25, but 
the measured values of Aqc and Bqc are not consistent 
with the theoretical predictions Moreover, we ob- 

serve such a linearity even in the 2D classical Heisenberg 
model where there should be no such crossover. Our data 
are not consistent with the one-loop order equations for 
the crossover regime 0H[§ either. Thus our results seem 
incompatible with the scenario of a crossover in £. 

The peak value of the static structure factor: Theory 
Till predicts for low T 



Sq(T) 

em 



= A2ttAT 



( T 



1-C- 



T 
l-nps 



(9) 



with a universal constant A. M w 0.307 [18| is the 
ground state magnetization. However the experimental 
data suggested ^,|| that Sq(T)/^ 2 (T) is temperature in- 
dependent over the temperature range accessible in the 
experiment. 

We find that the data at £ > 8.4 fit with an acceptable 
X 2 / n df = 0.9 (see Fig. |). The estimated value of 
C = —0.6(1) suggests that the effect of the correction 
term is less than one percent for T < 0.02, which may 
be regarded as the asymptotic regime for Eq. (g). For 
the universal prefactor we get the estimate A w 4.0. A 
series expansion study jL2| got A^, 2 w 3.2 for spin S = 
1/2 and A^ rj 6.6 for spin S = oo. Our result clearly 




FIG. 3. Ratio of the staggered structure factor and the 
square of the correlation length Sq/£ 2 versus T. The soild 
line is a fit to the theoretical low-T prediction Eq. ^ . 

shows, that as conjectured in Ref. |12]| , these values do 
not agree because the models are not yet in the low-T 
scaling regime. As the fits are unstable in the sense that 
A increases as we leave out data at higher T in the fit we 
view our estimate as a lower bound. 

We wish to note that for £ > 2.37(1) our data also 
fits to Sq(T)/£ 2 (T) - T N with N ~ 1.85(3) with an 
acceptable value of x 2 /Ndf — 0.75. 

Our data are definitely in agreement with the theory 
but not with the analysis from experiment Q . The com- 
parison with experimental data will be discussed in more 
detail in a forthcoming publication |15j . 

The staggered susceptibility: In the classical high-T and 
in the low-T renormalized classical regime the staggered 
susceptibility is expected to be related to the staggered 
structure factor as 



'Q 



Txst 



(10) 



This prediction agrees with QMC measurements by 
Sandvik et al. p0| . Our data, shown in Fig. || also con- 
firm it, with even higher precision. In the quantum crit- 
ical regime this ratio is expected to be also constant but 
with a value of MS3 



Txsi 



1.09. 



(11) 



Instead of an expected plateau with this value the data 
however show a peak with a maximum around T ~ 0.8 
p0[ . While the peak value is close to the predicted quan- 
tum critical regime the temperature is outside the range 



3 



0.95 




renormalized classical 



P= 1/T 



FIG. 4. Ratio of the staggered structure factor and suscep- 
tibility T ^ Q T versus = J/T. This ratio should be 1 (dashed 
line) both in the renormalized classical (T — + 0) and high tem- 
perature regimes (T — > oo). In the quantum critical regime it 
is expected to be 1.09. 

where quantum critical behavior is observed in the uni- 
form susceptibility. Thus we hesitate to ascribe this max- 
imum to quantum critical behavior. 
Conclusion and discussion: We believe that we have pre- 
sented for the first time a set of thermodynamic data 
which displays the asymptotic behavior predicted by chi- 
ral perturbation theory. The values of the spin-stiffness 
constant and spin wave velocity agree well with those 
given by the spin wave theory of the QHA. This agree- 
ment is truly remarkable since the chiral perturbation 
theory is on the NLcxM whereas the spin wave theory is 
on the QHA. Our results strongly confirm the validity of 
the map from the QHA to the NLcM in describing the 
long distance behavior of the former. Our estimates of 
p s and c from the analysis of \ u may be compared with 
previous estimates based on the size dependence formula 
near the ground state, i.e., p s = 0.185(2) and c = 1.68(1) 
jL3[. They are extracted, however, by three parameter 
fits, which generally involve rather large uncertainties in 
estimates. 

The fact that the asymptote of \u manifests itself only 
at such low temperatures may account for previous puz- 
zles raised in studies of the correlation length. The de- 
viation of £ from Eq. ([?]) is reduced to approximately 10 
percent at j3 ~ 4.9 from the 20 percent deviation seen 
previously at j3 = 2.5. In fact this is similar to the 2D 
classical Heisenberg model as demonstrated by recent nu- 
merical studies pj|)p2"[ . Our data of Sq and Xu also in- 
dicate that one needs to probe data with much lower T 
for the correction term to be safely ignored. 

Another puzzle raised in experimental ||,|| and nu- 



merical |12|] studies was: why is the agreement of the 
measurements worse for larger spins S? Equation (0) 
depends on the value of S only implicitly through the 
(S-dependence of p s and c. The reason for the larger 
discrepancies at larger S is implicit in the assumptions 
in the theory: it is valid only for « 1 and 

tt~ — << tt~ ~ 1/S. The latter condition restricts 
the validity of the theory to larger and larger correlation 
lengths as the spin S is increased. 

We also clarify the issue of the existence of the 
crossover. As suggested in Ref. j(| the uniform suscep- 
tibility, where all logarithmic corrections in the theory 
cancel, shows a clear crossover from a quantum critical 
to the renormalized classical regime around T ~ J/3. 
In all other quantities however no good evidence for a 
crossover can be observed. 

After completing our simulations, we became aware of 
a recent related preprint |23[ that addresses similar sub- 
jects. The authors claim to extract thermodynamic £ up 
to (3 = 12, with corresponding £ of order of 10 5 , based on 
recently developed finite size scaling extrapolation meth- 
ods plfl . This finite size scaling method requires the ex- 
istence of a universal scaling function that has no explicit 
T dependence. The existence of this function can be ex- 
plicitly checked numerically. However, contrary to the 
corresponding classical case pl| , we observed, in a high- 
precision study, that in this quantum spin system the 
scaling function becomes approximately T independent 
only when T is very small. Detailed accounts of this is- 
sue will appear in a separate paper 0] . 

The simulations were performed on the 1024-node Hi- 
tachi SR2201 massively parallel computer of the com- 
puter center of the University of Tokyo and used about 
250 000 hours of CPU time. M.T. was supported by the 
Japanese Society for the Promotion of Science. 



[1] Y. Endoh et al, Phys. Rev. B 37, 7443 (1987); K. Yamada 

et al, Phys. Rev. B 40, 4557 (1989); 
[2] M. Greven, R. J. Birgeneau, Y. Endoh, M. A. Kastner, 

M. Matsuda and G. Shirane, Z. Phys. B 96, 465 (1995). 
[3] K. Nakajima, K. Yamada, S. Hosoya, Y. Endoh, M. 

Greven and R.J. Birgenau, Z. Phys. B 96, 479 (1995) 
[4] S. Chakravarty, B. I. Halperin, and D. R. Nelson, Phys. 

Rev. B 39, 2344 (1989) 
[5] P. Hasenfratz and F. Niedermayer, Z. Phys. B 92, 91 

(1993) 

[6] A. V. Chubukov, S. Sachdev, and J. Ye, Phys. Rev. B 
49, 11919 (1994) 

[7] F. D. M. Haldan e, Phys. Rev. Lett . 61, 1029 (1988). 

[8] M. Troyer et al, |cond-mat/9702077| , J. Phys. Soc. Jpn., 
in press; M. Troyer and M. Imada in Computer Simula- 
tions in Condensed Matter Physics X, ed. DP. Landau 



4 



TABLE I. A selection of thermodynamic data for the 2D QHA in the range 0.25 < (3 < 4.95. 
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0.25 


1.75 


2.25 


2.75 


3.25 


3.75 


4.25 


4.75 


4.95 


e 


0.289(2) 


2.37(1) 


4.46(1) 


8.38(1) 


15.7(1) 


29.0(1) 


52.8(2) 


95.7(3) 


120.5(4) 


Xst 


0.0811(2) 


4.185(7) 
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1082(3) 
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4980(17) 


Xu 
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0.0735(1) 


0.0656(1) 


0.0604(1) 


0.0572(1) 
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-0.6432(1) 


-0.6559(1) 


-0.6618(1) 
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-0.6663(1) 


-0.6676(1) 




0.329(1) 


2.63(1) 


5.96(1) 


14.40(2) 


36.91(5) 


95.7(2) 


254.8(5) 


682(2) 


1007(5) 
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